Threshold features in transport through a ID constriction. 
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Suppression of electron current A7 through a ID channel of length L connecting two Fermi liquid 
reservoirs is studied taking into account the Umklapp electron-electron interaction induced by a 
periodic potential. This interaction causes Hubbard gaps Eh for L — > oo. In the perturbative 
regime where Eh -C v c /L (v c : charge velocity), and for small deviations Sn of the electron density 
from its commensurate values —AI/V can diverge with some exponent as voltage or temperature 
V,T decreases above E c = max(v c /L,v c Sn), while it goes to zero below E c . This results in a 
nonmonotonous behavior of the conductance. 
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Recent developements in the nano-fabrication tech- 
nique have made the ID interacting electron systems an 
experimental reality, and its quantum transport prop- 
erties have been the subject of extensive studies both 
experimentally and theoretically [3-15]. In realis- 
tic experimental set-ups, the quantum wire is attached 
to two-dimensional regions called reseivoirs or leads. To 
describe ID transport phenomena in this configuration 
a model was recently formulated of an inhomogeneous 
Tomonaga-Luttinger liquid (ITTL) |j|§]. It recovers 
the conductance G = 2e 2 /h observed experimentally 
even in the presence of the electron-electron interaction 
in the wire [ ^| , although the previous calculations on an 
infinitely long wire predicted the renormalized con- 
ductance. ( see [0 for the futher development). The same 
puzzle holds in the Integer Quantum Hall transport, c.f. 
H . Calculation of the conductance suppressed by a weak 
random impurity potential in this model Jl0| had agreed 
with both the previous theoretical prediction |7| and an 
experiment Q| (some amendment to that result will be 
pointed out below). Furthermore Tarucha et al. |l(| suc- 
ceeded to introduce the ID periodic potential with a pe- 
riodicity of order 40nm into the wire 2pm in length and 
50nm wide. This induces the Umklapp scatterings. The 
electron density n can be continuously controlled by the 
gate voltage, and one can satisfy the half-filling condi- 
tion within an accessible value of n. If this condition is 
satisfied the system will becomes a ID (doped) Mott in- 
sulator with the Hubbard gap Eh for the infinite length 
wire. Then it will offer an idealistic system to study the 
quantum transport in Mott and doped Mott insulators 
in ID. 

Inspired by these works, we study in this letter the- 
oretically the I — V characteristic of the wire of length 
L connected to leads taking into account the effects of 
the Umklapp electron-electron interactions, which can be 
directly compared with the experiments. We calculate 
the suppression of the current AI perturbatively in the 
Umklapp scatterings, and the results are summarized in 
Figs. 1 and 2 for a threshold structure near half filling. 



There are two energy scales, i.e., the finite size energy 
Tl = v c /L [v c : charge velocity) and E t hr ~ v c Sn [5n : 
the deviation of the electron density from its value at the 
filling v equal to 1/2) measuring the incommensurability. 
For T,V > E c — max(v c / 'L, Ethr)> the suppression of the 
current —AI/V diverges as max(T, V) i9 ~ 3 if the short 
range interaction constant g for the forward scattering is 
less than 3/4. For T, V < E c , on the other hand, the sup- 
pression —AI/V goes to zero as T, V — > 0. Then we pre- 
dict the nonmonotonous temperature and/or voltage de- 
pendence of /, which is the clear signiture of the Umklapp 
scattering effect. For small values of g, expected when 
the screening length £ c of the interaction is much larger 
than the width of the channel d and g oc \/ y/m S, c /d pA\ , 
we predict a few more threshold singularities. These fea- 
tures could be observed experimentally by changing the 
gate voltage, bias voltage, and temperature. 

Our model can be derived following pj from a 1 chan- 
nel electron Hamiltonian 



n = J dH£^)(- 
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E F )ip a (x) 



+ip(x)p 2 {x) + [V imp (x) + V perio d{x)]p(x)} (1) 

with the periodic potential V per i d(x) ( period a) as- 
sumed to be weak enough to justify the perturbative 
consideration of the Umklapp backscatterings. The 
Fermi momentum kp and the Fermi energy Ep is de- 
termined by the filling factor v as v — kpa/n and 
Ep w vpk F . In Eq. (||) the function ip(x) — 9(x)9(L — x) 
switches on the electron-electron interaction inside the 
wire confined in < x < L. Contribution of the 
random impurity potential Vjmv(x)p(x) to the conduc- 
tance has been considered in 0,|lO[, some results of 
which we will use below. Following Haldane's general- 
ized bosonization procedure [fLsf to account for the non- 
linear dispersion one has to write the fermionic fields as 
ip a {x) = v / k F /(2TT)Y / exp{i(n + l)(k F x + <j> a {x)/2) + 
i9,y(x)/2} and the electron density fluctuations as p(x) = 
J2pcr(x), p a (x) = (d x 4>< J (x))/{2ir)Y J exp{in{k F x + 
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(j} a (x)/2)} where summation runs over even n and <p a , 6 a 
are mutually conjugated bosonic fields [(f) <7 (x),0 <7 (y)} = 
i2itsgn(x — y). 

After substitution of these expressions into (|l|) and 
introduction of the charge and spin bosonic fields as 
<\>c.s = {4*] i <Al)/\/2 the Hamiltonian takes its bose-form 
TL = TLo + Tibs ■ Here the free electron movement modi- 
fied by the forward scattering interaction is described by 
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with g c (x) = g for x G [0, L] ( g is less than 1 for the repul- 
sive interaction and it will be assumed below ), g c (x) = 1, 
otherwise and v c (x) — vp/g c (x). The constants in the 
spin channel g s = l,v s = vp are fixed by SU(2) symme- 
try. Keeping only the most slowly decaying terms among 
others with the same transfered momentum one could 
write the backscattering interaction as 



nj E F 
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VF Jo 



dx 



E 

even m>0 



U m cos(2k m Fmx 
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Difference in the transfered momentums 2k m p is brought 
up by the periodical potential with the period a: k m p — 
kp — irl/ (ma)) > 0, where I is an integer chosen to mini- 
mize k m p . We have omitted the first terms in both sums 
for they could not affect the current in the lowest pertur- 
bative order: the m = term of the first sum can contain 
only the spin field and the m = 1 term of the second sum 
is assumed to transfer large momentum. The dimension- 
less coefficients U m originating from if(x) will be assumed 
to be small enough to justify perturbative calculation of 



the current. 

Change of current due to the backscattering is 
given by : Af = -i/(2y/2w) J dc [d x 9 c (x),H] = 
\[2 J dx(5/Scj) c (x))Hbs- At finite voltage V applied sym- 
metrically to neglect the momentum transfer variation, 
the average of AI decomposes into sum of the different 
backscattering mechanism contributions < AI m > in the 
lowest order. The even m terms involving only <fi c field 
are equal to 



m 

< Af TO >= - — 



• U m Ep % 
' vf ' 



dt 



dr\dx2 < e 



im(p c (x\ f t)/v2 f> — im(p c (x2 ,0)/v2 ^ [^im(2k m F {x\ — X2 )+Vt) 



h.c 



(4) 



The current operator AI m has a high energy scaling di- 
mension m 2 g/2 and a free electron (g — 1) behavior at 
low energy. We will see below that the integral (|]) scales 
at low energy with (m 2 — 1) exponent and with (m 2 g — 2) 
exponent at high energy. The most singular behavior is 
due to Umklapp backscattering at m = 2 with the thresh- 
old voltage V = E t hr = 2k2pv c going to zero at the half 
filling. Less singular correction with m — 4 could become 
relevant at the one and three quaters fillings and so on. 
Expressions for the odd m terms include additionally a 
spin field correlator < e i<t>s{xut)/V2 e ~i<t> s (x2.,o)/V2 > under 



the integrals in (0). The high energy dimension of M m 
in this case is m g/2 + 1/2. The most singular behavior 
occurs to the m = 3 term at the one and two thirds fill- 
ings. It has two threshold energies E t hrc.s = 2k^pv c ^ s for 
v c ^ v s . 

charge field exponents e 1 ^"^'^ , 
is specified by Ho, could be 
correlators of the uniform TL 
K(x,t,g,v c ) (K(x,t,g,v 



Correlator of the 
evolution of which 
compiled from the 
liquid K(x,t) = 



(W/3) 2 7( 

lowing way 



sinYi 9 (-K{x/v ±{t - ia))/P))) in the fol- 



< e i<t>c(x,t) e -iMv,o) >= K(x-y,t) Yl ( 



K(2nL,0) 



l^K(2nL±\x-y\,t)' 



^ ( K(2(nL + x),0)K(2(nL + y),0) - r - 
K 2 (2nL + x + y,t) ' 



n=0 



72 fr ( K(2(nL - x), 0)K(2(nL - y), 0) ^-^/2 
K 2 (2nL~x-y,t) ' 



71=1 



(5) 



Here (3 is inverse temperature 1/T and a — I/Ef is 
the ultraviolet cut-off. This complicated form comes 
about through a multiple scattering at the points of joint 



x = 0, L. As a result of the scattering the correlator 
< cj) c (x,t)(l) c (y,0) > becomes an infinite sum of the uni- 
form correlators taken along the different paths connect- 



2 



ing points x and y and undergoing reflections from the 
boundaries at x = 0, L. Each reflection brings addi- 
tional factor r = (1 — g)/(l + g). The similar correlator 

< el 4> s (x,t) e - lt p s (yfi) > for gpm ' fidd ig K ( x _ y^\,V a ). 

Below we analyze the current corrections (|J) for high 
(T > l/t L = T L ) and low (T < T L ) temperatures, re- 
spectively. 

I. High temperatures T > Tr, - The uniform correlator 
K{x, t) goes down exponentially if distance between the 



points | x | exceeds the inverse temperature. Therefore 
only paths with length less than (3 contribute to the corre- 
lator (Jfy. This means that the high temperature form of 
the correlator (||) reduces to the first multiplier K(x—y, t) 
up to a factor (1 + 0{exp{-t L / (3)). Neglecting 0(T L /T) 
quantity we can extend integration over x\ — X2 in (^) 
from — oo to +oo. Then calculation of the m = 2 contri- 
bution reduces to finding Fourier transformation F 2g (q, e) 
of the correlator K 2 (x,t, g,v c ): 



< Af 2 >= - 



1 [U 2 E 



t L TF 2g (2E thr ,±2V) = -2( 



T(2g)g > Tr, Ep 



V±E thr s 
2-kT ' 
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One can easily see its behavior making use of the following asymtpotics: 

2 P 2 



< Af 2 >oc 



" V 1 Tl { ({V + E thr )T/E F fB-\ V « £ t/ir » T 



2^2 



(7) 
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FIG. 1. Schematic voltage dependence of the high temper- 
ature current corrections produced by the m = 2 Umklapp in- 
teraction A/2- Solid lines, E t hr = 0; dashed lines, E t hr T. 



These asymptotics show that the threshold singu- 
larity in the current voltage dependence diverges as 
{(V - EthrVT) 23 - 1 if g < 1/2 (Fig.|). It becomes 
stronger in the differential conductance dependence. At 
g > 1/2 the differential conductance correction <G 2 be- 
haves as -(V/E F ) 4 s- 3 , and saturates at -(T/E F ) 4 3- 3 
below T, if E t hr < T; otherwise, the correction shows 
divergence —((V — E t hr) /T) 29 ~ 2 smeared over T scale 
near the threshold and becomes suppressed exponen- 
tially as —exp(~E t hr/T) below it. Generalization to the 
other even m current corrections needs just changing: 
4g — > m 2 g , Ak 2F ± 2V — > m(2k mF ± V"). The edge singu- 
larity is characterized by a half of the scaling dimension 
for AI m since only one chiral component of the field <f> c 
contributes. As to the odd m terms, the two treshold en- 
ergies Ethrc.s = 2k mF v c . s become distinguishable if their 
difference exceeds T. The leading high-temperature cur- 
rent correction reads as: 



< Af m >= — 



U, 



iE F 



2ng 



J^T J J dqfeFi {mEthrs ~ qv Sl ±mV - e)F m 2 g (qv c ,e) 



(8) 



Substituting zero temperature form of F a func- 
tion F a (q,w) = 8[sin(na)T(l - a)] 2 (a/2) 2a f[ ± (w ± 
q) a ~ 1 8(w ± q) in Q one can gather that the current cor- 
rection behaves as —(V/E F ) m 9 ~ 1 at large voltage V > 
Ethr c,s, has a leading singularity -((V - E thrc )/T) m2g / 2 
smeared over T scale near the first threshold and — ((V — 
Ethr s)/T) m 9 ~ 1 / 2 near the second threshold (we assume 
v c > v s ). Below the lowest threshold it becomes expo- 
nentially suppressed. These singularities result in the di- 
vergences of the differential conductance or higher deriva- 



tives of the current in voltage. The threshold behavior of 
the m = 3 term of the differential conductance correction 
is divergent at Ethr c if 5 < 2/9 and at E t h r s if 5 < 1/6. 

2. Low temperatures T,V << Tl - With lowering tem- 
perature we should expect that above current correction 
dependences will be modulated by a ttTt, quasiperiodical 
interference structure |l4|,[l5| and also a new low energy 
scaling behavior of the current correction operators will 
appear at V/T < Tj, . The dominant contribution to the 
integral of (0) comes from long times t ^S> t F . One can 
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neglect the spacious dependence compared with large t tions n < n* — /3/(2i£7r) only to come to the long time 
in (Q) and keep the multipliers with the number of reflec- asymptotics: 



< e 



■Mx,t) e -t4, c ( v ,Q) >= e 7(rz,/T)/^.\2s/ (irt L //3) 2 > 

ti, sinh(w(t — ia) / f3)sinh(ir(—t + ia) / (3) ' 



. y/xy{L- x)(L-y) . 2r 
L 2 > 



(9) 



where z{T L /T) = r ^ (tL ^ and j(T L /T) approach the 
constant 7(00) on the order of 1 as \n(T L /T)z(T L /T). 
Our asymptotic analysis following in essential Maslov's 
paper |l(fl shows that the low energy exponents approach 
their free electron values as exp[Ti, \xir / (T-k)\. The effect 
accounts for prolongation of the paths due to the finite 
reflection. In particular, it determines the coefficient c(g) 



of the T 2 corrections to the non-universal zero tempera- 
ture value of the conductance variation due to impurities: 
AG lmp oc -{L/I^T/Ef)^ 1 ^ - c(g){T/T L ) 2 ) in a uni- 
versal way |p7| . After substitution of (||) into Eq. (Q), 
the current suppression produced by the even m terms 
of the interaction becomes equal to: 



l2 m 2 (l- Z ) e m 2 7 /2 /jj 



r(m 2 (l - Z)) 



!1 



R m 2 g (2mk mF L)T L 



ttT 



E f 



2 (9-1) 
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where function f a (x) = sinh(a;) Y[± r(a ± ix/n) characterizes the V — T cross over. It approaches T 2 (a)x(l — 
(lnr(a))*(a;/7r) 2 ) at x <§; 1 and ■K(x/ir) 2a ~ 1 at x ^> 1. Function R specifies the kp — 1/L cross over as: 

R2 a( x ) = ji+trg J i/2+2r 9 (^/2) ^ T 2 (l + 2rg) I ll^( x /2 T - wIg)x- 2 - g ^\ I > 1 




FIG. 2. Schematic temperature dependence of the conduc- 
tance corrections produced by the m — 2 Umklapp interaction 
AG2 at g < 3/4 (linesl) and at g > 3/4 (lines2) and by the 
random impurity potential AGi mp [|7jJlC|] (lines 3) . Solid lines, 
E t hr = 0; dashed lines, E thr S> T L . The dot lines 1 and 2 (1' 
and 2') are the full conductance correction at Ethr = (finite 

Ethr)- 

The odd m corrections will meet Eq.(|l(]) after substi- 
tution m 2 + 1 instead of m 2 . Combining above results 
we can outline a temperature dependence of the conduc- 
tance correction produced by the m = 2 Umklapp in- 
teraction (Fig.^). Its magnitude increases/decreases fol- 
lowing (E F /T L )(T/E F ) 4 9- 3 as T going down above T L 
and follows (T/T L ) 2 (T L /E F ) 4 3~\ if E thr < T L ; oth- 
erwise, the correction starts to decrease exponentially 
exp(—E t hr/T) below E t h r and keeps on decreasing like 



(T/T L ) 2 (T L /E F )^-\T L /E thr ) 2 +^sin 2 (2k 2F L - nrg) 
below Tl. The T(> Tl) dependence is similar to that 
of the conductivity of infinite wire found by Giamarchi 
fL2fl . Similar dependence with T replaced by V could be 
predicted for the zero temperature differential conduc- 
tance dG 2 (V) at V < T L . 

In summary under perturbative condition we have de- 
scribed a hierarchy of the threshold features produced 
by the Umklapp backscatterings at the rational values 
of the occupation number inside the ID channel con- 
necting two Fermi liquid reservoirs. In the differential 
conductance (its derivative) vs. volatage, the threshold 
structure is an asymmetric peak of width max{T,Ti,) 
located at treshold, E t hr- In the conductance vs. tem- 
perature, we predicted a maximum below E t hr due to 
crossover from the Umklapp backscattering to the impu- 
rity suppression and an asymmetric minimum at E t hr if 
the interaction is strong enough. 
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